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ABSTRACT 

We present the near-infrared images and spectra of four silhouette disks in the 
Orion Nebula Cluster (M42) and M43 using the Subaru Adaptive Optics system. 
While d053-717 and dl41-1952 show no water ice feature at 3.1 fim, a moderately 
deep (xjce ~ 0.7) water ice absorption is detected toward dl32-1832 and d216- 
0939. Taking into account the water ice so far detected in the silhouette disks, the 
critical inclination angle to produce a water ice absorption feature is confirmed to 
be 65-75°. As for d216-0939, the crystallized water ice profile is exactly the same 
as in the previous observations taken 3.63 years ago. If the water ice material is 
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located at 30 AU, then the observations suggest it is uniform at a scale of about 
3.5 AU. 

Subject headings: dust, extinction - evolution - infrared: ISM - infrared: plan- 
etary systems - protoplanetary disks - stars: individual (d053-717, dl32-1832, 
dl41-1952, d216-0939) 



1. Introduction 



A silhouette disk is one of the most suitable targets to investigate the primordial com- 
position of a planetary disk by using its central star as a background light source. The Orion 
Nebula Cluster (ONC; M42) and its neighbor H II region M43 are the nearest high-mass 
star-forming regions, exhibitin g a number of silhouette disks. Since the first discovery by 
McCaughrean &: O'Delll (Il996[). the number of the silh ouette disks has significantly increased 



flBallv et al 



2000; 



Smith et al.ll2005t iRicci et al.ll2008l ) and currently the number of the pure 



silhouette disks is 
features. 



■^30, which enables us to conduct a systematic study for disk absorption 



Water ice in protoplanetary disks is particularly important, because it is thought to 
evolve to icy bodies such as comets in our solar system. In theory, the distribution of the 
water ice in the protoplanetary disks has bee n extensively investigated as a key parameter 
for the planet formation (e.g., iHayashil Il98l[ ). which predicts the inner boundary around 



2 AU of the water ice presence ("snow line") where the temperature is 170K for the ice 
condensation. Although the mid-plane of the disks is of the greatest interest for the planet 
formation, i t is difficult to obseryatiori ally probe into the mid-plane of the disks due to high 
extinction. iPontoppidan et al.1 ( l2005l ) constructed the model of the water ice absorption 



in the circumstellar disk through the water ice detection for the edge-on disk object CRBR 
2422-3423 in p Ophiuchus dark cloud, which suggests a critical inclination angle of the disk for 
strong water ice detection to be around 69°. Since silhouette disks show a clear morphology 
of the circumstellar disks, the inclination angle of the disks can be more accurately estimated. 
In fact, the silhouette disks located in M42 and M43 show a wide range of the disk inclination 
angle, and those are very suitable for examining the critical inclination angle in a larger 
sample. The silhouette disks in M42 and M43 are exposed to the intense UV radiation 
from O and B stars, and the critical inclination angle may be increased to a larger value 
by dissociation and desorption process of the outer layer of the disks. Since star formation 
occurs more commonly in high-mass star-forming regions, it is important to understand the 
water ice formation under such environments. 
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Strong water ice absorption at 3 /xm band is very useful to investigate characteristics 
of the water ice materials. In particular, it is well known that the profile and peak of the 
water ice absorption in the wavelengths of 2.9-3.3 fim provides the information of its grain 
size and crystallinity. While larger water ice grains bring its absorption peak t oward longer 



wave length, more crystallinity causes a sharper peak at longer wavelength ( ISmith et al. 



19891 ). The water ice feature so far detected at 3/im band toward young stellar objects is 
usually explained by the combination of amorphous and crystallized water ice absorption 
with a grain size of 0.1-0.5 /im, whose pe aks are located at 3.0-3.05 /zm and ~ 3.1 /zm, 
respectively ( iDartois fc d'Hendecourtll200ll ). 



Terada fc Tokunagal (120121 . hereafter Paper I) carried out L-band spectroscopy of five 



bright young stellar objects with silhouette disks in M42 and M43 and detected a deep 
water ice absorption (rjce=0.67) toward d216-0939 in M43, which is associated with a large 
silhouette disk (~1000 AU in diameter). In the profile of the water ice absorption, the 
3.2 Atm enhancem ent that is characteristics of crystallized water ice with large grains (~1 fim; 
Leger et al.l 119791 ) was seen, and it su ggested grain growth a nd thermal processing in the 
disk. Since the nearby star JW 671 1 Jones fc Walkerl Il988h shows no water ice feature, 
the water ice detected for d2 16-0939 is concluded to reside in the protoplanetary disk of 
d216-0939. Whereas d l21-1925 also shows the water ice absorption, its nearby star JW 370 
( iJones &: Walkerlll988l ) exhibits the same depth and profile of the water ice absorption, which 
suggests that the water ice toward dl21-1925 is in the foreground. There is no evidence of 
the water ice feature detected in the remaining three objects: dl82-332, dl83-405, and d218- 
354. The low inclination angle of the disk is likely the primary cause for no water ice feature 
in the disk in the line of sight for dl21-1925, dl82-332, dl83-405, and d218-354. Another 
possible cause is low abundance of the water ice in the disk due to the strong UV radiation 
from O and B stars. Increasing the number of spectra will help to understand the prevalence 
of water ice in the silhouette disks. 

Adaptive optics (AO) with large ground-based telescopes gives us a unique opportunity 
to obtain the highest spatial resolution with diffraction limited images, whic h reaches to < 
O^.^l resolution at 3 urn. At the distance to M42 (414 pc; iMenten et al1l2007h and M43 (440 
pc; lO'Dell fc Henneyll2008l ). the direct study of the outer region of protoplanetary disks can 



be done with 40 AU resolution. Many circumstellar disks have been so far detected using 
AO, including t he disks with high inclination angle (>80°) which sometimes show an e dge- 
on morphology (IMonin fc Bouvieii l2000l : iJayawardhana et al.l 120021 : iPerrin et al.l |2006| ) . In 
addition, AO observation has another significant advantage to provide increased sensitivity 
for photometry of point sources and use of a narrower slit for spectroscopy. It is crucial 
to use AO to increase the reachable targets associated with the silhouette disks. For AO 
observations, a nearby and bright reference source to sense the wavefront of light from the 
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target region is required. For targets that do not have such nearby and bright natural stars, 
an artificial laser guide star (LGS) is available at some large telescopes. In the outskirt 
regions of M42 and M43, LGS plays an essential role to make the observation of silhouette 
disk possible. 

Using the Subaru Telescope, we present natural guide star (NGS) and LGS images and 
spectra of four sources with silhouette disks in the outskirt of M42 and M43 regions, na mely 
d053-717, dl32-1832, dl41-1952, and d216-0939 JSallv et all [2OO0I : ISmith et all boosh . In 
Section 2, our sample selection is explained. The observational procedures and applied data 
analysis are described in Section 3. In Section 4, the AO images and spectra for all the 
targets are shown, and the water ice optical depths are derived after the spectral continuum 
fit. The water ice feature in each disk is discussed in Section 5. The summary is described 
in Section 6. 



2. Samples 



2005 



Targets wer e selected from all the pure silhouette disk samples (ISmith et al 
Ricci et al.l 120081 ) with three criteria of (1) an inclination angle of >55°, (2) a distance of > 



200" from the O and B stars, and (3) a brightness of K' < 15.5. As for the first criterion 
about the inclination angle. Paper I found that the silhouette disks likely have the critical 
inclination angle of 65°-75° for the detection of the water ice in disks, and it is required to 
confirm this possible critical angle by observing more targets with inclination angle 65°-75°. 
The second criterion is to avoid the environmental effect of strong UV radiation from the 
O and B stars. In fact, the previous samples with no water ice feature in Paper I are all 
close to the Trapezium cluster, and it was difficult to distinguish the environmental cause 
and the geometrical cause for no water ice detection. The third criterion is for obtaining 
sufficient signal-to- noise ratio (S/N) for spectroscopy at L. Typically, L < 14 is required 
for detection of absorption feature with good S/N using an 8 m telescope. We selected the 
targets which met this criterion through a K' imaging survey (H. Terada et al., in prepara- 
tio n). In Figure [H the loca tion of our targets is superposed on the extinction map derived 
by IScandariato et al.l (120111 ). Table [1] summarizes the disk parameters of each target. 



3. Observations and Data Reduction 



All the observations were perfornied us i ng Infrared Camer a and Spectrograph (IRCS: 
Tokunaga et al.lll998l : iKobayashi et al.ll2000l : iTerada et al.l 120041 ) combined with the Subaru 
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Table 1. Disk Parameters 









Inclination 


Exciting 




Object 


Mass*^ 


Diameter 


Angle 


Source 


Distance 




(xlO-2 Mq) 


(") 


n 




(") 


d053-717 


0.52±0.07 


0.9^ 


65-85^ 


0^ Ori C 


286.32 


dl32-1832 


0.80±0.08 


1.5" 


75" 


6^ Ori C 


294.19 


dl41-1952 


1.50±0.06 


0.7^ 


55-60'^ 


6^ Ori C 


214.18 


d216-0939 


4.53±0.06 


2.6^ 


75-801^ 


NU Ori 


400.49 



Mann fc Williamsl f l2010l ). 



Smith et al. 


(2005) 


'Ballv et al. 


2000). 



Note. — Distance is measured from the listed 0,B stars responsible for 
excitation of the H II region. 
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Fig. 1. — Location of our samples. Asterisks represent the exciting O and B stars, namely 
NU Ori for M43 and 6^ Ori C for M42. Open squares are the previous samples in Paper I, and 
open circles are corresponding to the samples in this paper. d216-0939 is observed in both 
runs. The sampl es are superposed o r i the visual extinction (Ay) map for the ONC members 
from Figure 6 of IScandariato et al.l (120111 ). Gray and black contours are at the Av=3 and 
6 mag, respectively. d053-717 is found to be located in the region of the least extinction {Ay 
< 3) among the four targets. dl41-1952 and d216-0939 have moderate extinction {Ay ~ 3), 
and there is the largest extinction around the region of dl32-1832 {Ay ~ 10). 
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Adaptive Optics system (A0188: iHayano et al.ll2010l ). IRCS provides functions of imaging 
and spectroscopy with grisms and a cross-dispersed echelle in the wavelength range of 0.9- 
5.6 fim, and it is optimized for the high angular resolution images delivered by the AO. The 
observing log is summarized in Table [21 which includes information about the airmass of each 
observation and standard stars. The NGS mode was used for observations of d216-0939, since 
a sufficiently bright star (i?~10 mag) for the wavefront sensing is located at 23'.'42 north away 
from the object. For the o ther targets: d053-7 17, dl32-1832, and dl41-1952, the observations 
utilized the A0188 LGS ( iHayano et al.ll2010l ) for higher order wavefront correction together 
with a tip-tilt and defocus correction using an i?~14 mag star. As for d053-717 and dl41- 
1952, the targets themselves are bright enough (i?=13.7 mag and i?=13.6 mag, respectively) 
to be used as a reference star for the lower order of wavefront correction (tip-tilt guide star, 
TTGS). In the case of dl32-1832 observations, the bright star (i?=14.8 mag) located at 49" 
southwest away from the object was used as TTGS. 

Imaging was conducted by using a nine-point box-shaped dithering with an offset sep- 
aration of 5'.'0. Spectroscopy was performed using A-BB-A dithering sequence along the slit 
with an offset separation of 1'.'5. The slit width was 0'.'225, which corresponds to a spectral 
resolution (A/AA) of 850 at K and 510 at L. Sky flat was used for L', and the fiat fielding for 
the J, H, and K' images and K and L spectra is done with an integration sphere illuminated 
by a halogen lamp. 

Standard analysis procedure is applied using the IRAF packages for all the data reduc- 
tion of sky subtraction, fiat fielding, shift-and-add and combining of the images and spectra, 
telluric correction, and w' avelength calibration. FS 13 (ii''=10.126) and FS 154 (ii''=11.050; 
Casali fc Hawarderu 119921 ) were used as K' photometric standard stars. For L' photometric 



calibration, observations of HD 40335 (L'=6.441) and HD 1160 (L'=7.055; iLeggett et al. 



20031 ) were conducted. For correction of the telluric absorption, HR 2075 (AOIV), HR 2328 
(AOVn), and HR 1826 (AOVn) were observed with an assumption that those stars have a 
blackbody continuum with a temperature of 9480 K. The hydrogen absorption o f these spec- 



trosco pic standard stars is appropriately removed using the method developed by lVacca et al. 



( 120031 ). The Br 7 emission at 2.166 /xm is likely the result of poor subtraction of the nebular 
emission and is not intrinsic to the source. The airmass mismatch between the object and 
the standard observations was very small as shown in Table |5J Wavelength calibration was 
done using the telluric absorption lines in the spectrum. 
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Results 



4.1. 1.2—4.2 yum Images with Adaptive Optics 

Aperture photometry is apphed for all the images with a radius of l'/5, and the photo- 
metric results are summarized in Table [31 These photometric values include the scattered 
light component around the central star. 

Images of the targets are presented in Figure [21 which show good image quality with 
FWHM of <0'.'2. Table[llshows the achieved spatial resolution and Strehl ratio for the images. 
The Strehl ratio is calculated from the peak intensity normalized by that of point-spread- 
function (PSF) without any turbulence. To investigate possible extended features around 
the sources, subtracted images by their nearby bright stars are presented in the bottom 
panels of Figure [21 The reference stars for the PSF are chosen at a vicinity of the objects. 
JW 283, 2MASS J05351361-0518329, JW 413, and JW 671 are located at 10'.'79 northeast, 
5'.'68 east, 7'.'04 west, and 4'.'38 southwest away from each object, respectively. Whereas it 
is well known that the PFS obtained with an AO system is affected by the anisoplanatism, 
the Strehl ratio degradation for the sep aration of 10^^ is only up to ~ 0.1 by n ieasuring with 



the A 0188 system in both NGS mode ( IMinowa et al.|[2010[ ) and LGS mode ( IMinowa et al. 



2012[ ). The anisoplanatism for the PSF reference stars with an offset of < 10'.'79 does not 
affect the following analysis. In the subtracted images, extended features are seen around 
K' image of dl32-1832 and J, H , K', and L' images of d216-0939. As for the images of 
dl41-1952, the reference star is fo und to be extended, w hich is known to be a young stellar 
object with a fluorescent Fe line ( iTsujimoto et al.[[2005[ ). As a result, a negative extension 
appears at both K' and L'. 

The extended feature at K' of dl32-1832 is elongated along the position angle of the 
silhouette disk. In geometry of nearly edge-on disks, the dark lane is usually seen between 
two sides of scattered hght from the central star. In fact, the Hubble Spac e Teles cope (HST) 
Advanced Camera for Surveys (ACS) image at Jc, shown in [Ricci et al.[ ( 12008[ ). marginally 
exhibits such morphology. It is probable that the spatial resolution of 0'.'2 is not enough to 
resolve the possible dark lane of this object in the infrared. Also, no outflow signature of 
molecular hydrogen emission which is seen in our spectrum of dl32-1832 is exhibited in the 
subtracted image at K'. The object d216-0939 is known to have a la rge reflectioii nebul a 
extended by ~3" and was discovered in the optical using HST ACS by [Smith et al.[ ( |2005[ ). 
The infrared extended feature detected at J, H, and K' resembles the morphology in the 
optical. Regarding the reflection nebula at L', it is detected on only the east side of the 
dark lane and the surface brightness is much fainter than that detected in the J, H, and K' 
images. 
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Table 2. Observing log 





Observing 




Exposure 


Average 


Standard 


AAirmass 


PSF 




Object 


Date 




Mode 


Time 


Airmass 


Star 


Obj.-Std. 


Reference Star 


Offset 




(UT) 






(s) 










( ) 


d053-717 


2011 Oct 


23 


imaging 


15 


1.115 


FS 13 


—0.040 


JW 283 


10.79 




2011 Oct 


23 


L' imaging 


50 


1.118 


HD 40335 


-0.036 








2011 Oct 


23 


K spectroscopy 


120 


1.113 


HR 2075 


-0.026 








2011 Oct 


23 


L spectroscopy 


480 


1.108 


HR 2075 


-0.015 






dl32-1832 


2011 Oct 


23 


K' imaging 


90 


1.164 


FS 13 


-1-0.008 


2MASS J05351361-0518329 


5.68 




2011 Oct 


23 


L' imaging 


100 


1.173 


HD 40335 


-1-0.020 








2011 Oct 


24 


K spectroscopy 


480 


1.175 


HR2328 


-1-0.047 








2011 Oct 


24 


L spectroscopy 


3240 


1.119 


HR 2328, 


-1-0.031 






dl41-1952 


2011 Oct 


23 


K' imaging 


3 


1.138 


FS 13 


-0.018 


JW 413 


7.04 




2011 Oct 


23 


L' imaging 


50 


1.140 


HD 40335 


-0.013 








2011 Oct 


23 


K spectroscopy 


60 


1.139 


HR 2075 


-0.001 








2011 Oct 


23 


L spectroscopy 


120 


1.127 


HR 2075 


-0.0001 






d2 16-0939 


2008 Sep 9 


J imaging 


240 


1.175 


FS 154 


-0.005 


JW 671 


4.38 




2008 Sep 9 


H imaging 


270 


1.161 


FS 154 


-0.024 








2008 Sep 9 


K imaging 


270 


1.194 


FS 154 


-1-0.007 








2008 Oct 9 


L' imaging 


560 


1.322 


HD 1160 


+0.070 








2009 Oct 2 


L spectroscopy 


1000 


1.104 


HR 1826 


-0.007 







Table 3. Total Magnitudes 



Object 


J 


H 


K' 


L' 


d053-717 






10.81 


11.12 


dl32-1832 






15.06 


13.62 


dl41-1952 






10.20 


9.53 


d216-0939 


14.43 


13.58 


12.35 


10.74 



Note. — The typical 1 a uncertainties 
are 0.02, 0.02, 0.02, and 0.1 mag for J, H, 
K\ and L', respectively. 
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(a) 



d053-717 



K' 



L' 




JW 283 




HSTACS F658N 



















Obj-Ref 






0-.2 




Fig. 2. — Adaptive optics images of the four targets with the silhouette disks. For each 
figure (a)-(d), the top panel shows the object image, and the PSF reference star is presented 
in the middle panel. Images in the top and middle panels are normalized by a peak flux and 
those are displayed in logarithmic scale with black for 100% and white for 0.2% of the peak 
flux. The bottom panel for J, H, K', and L' shows the subtracted image of the objects by 
the PSF reference. This figure is displayed in linear scale with black for 10% and white for 



10% of th e peak flux. 



or c omparison, the bottom-right panel shows the HST ACS Ho; 



image from lSmith et al.l (120051 ). All the images have the same size of 1'.'75 x 1'.'75. 
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(b) 



dl32-1832 



2MASS 
J05351361 
-0518329 



Obj-Ref 



Fig. 2. — Continued. 



Table 4. Achieved Spatial Resolution and Strehl Ratio 



Object 


Band 


FWHM 


Strehl 






(") 


Ratio 


d053-717 


K' 


0.092 


0.12 




V 


0.11 


0.44 


dl32-1832 


K' 


0.20 


0.05 




L' 


0.16 


0.30 


dl41-1952 


K' 


0.11 


0.13 




L' 


0.12 


0.48 


d216-0939 


J 


0.13 


0.035 




H 


0.11 


0.086 




K' 


0.081 


0.25 




L' 


0.12 


0.55 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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4.2. 1.9—4.2 yum Spectra with Adaptive Optics 

Figure [3] shows the spectra of the four targets: d053-717, dl32-1832, dl41-1952, and 
d216-0939. The level of the continuum at K and L is determined by estimating the flux 
in the slit taken from the images shown in Figure |2j The same procedure for the spectral 
model fitting as described in Paper I is a pplied for all the spectra using IRTF spectral library 



( ICushing et al.ll2005l : iRayner et al.ll2009l ). and the best-fit model spectra are overplotted with 
gray lines in the same panel of Figure [31 The best-fit parameters are shown together with 
the previously known spectral type and Ay in Table [51 It is noted that the derived Ay could 
be underestimated due to the scattered light around the central sources. 

A broad absorption feature around 3.1 /im is seen in the spectrum of dl32-1832 and d216- 
0939. On the other hand, there are no absorption signature at 3.1/im band for d053-717 
and dl41-1952. The center wavelengths, optical depths, and band-widths of this absorption 
feature for each target are derived from the spectra (Table [B|). Whereas S/N for dl32-1832 
is relatively low (~10 at the continuum) due to its faintness at L' (13.62 mag), the profile 
of the 3.1 /im absorption matches the water ice band. In contrast, the nearby star 2MASS 
J05351361-0518329 shows no water ice feature in the spectrum (Figure [1]) based on the best- 
fit model spectra with a spectral type of M7V, Av=6 mag, and a blackbody temperature of 
1000 K. 

For the case of d216-0939, the detection of the apparently comparable absorption at 
3.1/im band on different epoch (2006-01-15-UT) has been reported in Paper I, and this 
water ice band is found to be unchanged after a time interval of 3.63 years. Since d216-0939 
is spatially resolved in the L' image, spectra are extracted for the three locations (northeast 
(NE), center, and southwest (SW) of d216-0939) shown in the left panel of Figure [3 Using 
the PSF of the nearby star JW 670, contribution of the wing of the PSF at L' to NE and 
SW is estimated to be 12% and 13%, respectively. Thus, the spectra for NE and SW are 
affected by the contribution from the center location. The position angle of the slit was set 
to the direction from the object to the nearby star JW 671, and the slit is not aligned to the 
disk direction. Whereas the NE position covers the L' scattered light region effectively, the 
SW position is around the mid-plane area with some scattered light. The derived spectra for 
each position are shown in the left panel of Figure [51 The photometric points at K' and L' 
are estimated from the flux density within the slit as measured from the AO images. Model 
spectra are determined using the K' photometric points and normalized spectra at L, and 
those are shown with gray lines in Figure [51 Both the positions of NE and SW show bluer 
continuum, and the water ice absorption is deflnitely more shallow, and the optical depths, 
center wavelengths, and widths of these spatially resolved water ice features are listed in 
Table [3 
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d216-0939 1 




HSTACS 









r.o 


— 



ACS F658N Jan 2ofl 










Apr 200S 








1".0 




ACS F6S8N 












1".0 





2 2.5 3 3.5 4 

Wavelength [ jLtm ] 

Fig. 3. — Spectra of the four targets associated with the silhouette disks. Black lines show 
the spectra of the targets, and gray lines are the model spectra. The absolute level of the 
spectrum is determined from the flux, which falls in the slit. In the r ight panel, HST ACS 
Ha images are shown with a size of 7" x 7" from ISmith et al.l (120051 ) . Only th e dl32-1832 
i mage in 2005 April shown as one of the inserts is the online data presented in iRicci et al. 
J2008h . 
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Table 5. Spectral Type, Visual Extinction, and Blackbody Temperature 



Hillenbrand (1997) Derived Parameter 





Spectral 


Spectral 




Blackbody 


Object 


Type Ay 


Type 


Av 


iemperature 




(mag) 




(mag) 


(K) 


d053-717 


K5-6 0.43 


M2±l 


1.5±0.20 


1100 


dl32-1832 




K5±2 


1.5±0.30 


1050 


dl41-1952 


Late G 2.55 


M2.5±l 


2.55±0.15 


1100 


d216-0939 


K5 0.72 


K5±l 


0.72±0.18 


1190 


Table 6. Optical Depth of Water Ice Absorption 


at 3.1 /xm 




Object Peak 




Az/ 






(/xm) 




(cm^^) 






d053-717 


<0.065 








dl32-1832 3.08 


0.70±0.05 


445 






dl41-1952 ■ ■ • 


<0.0052 








d216-0939 3.12 


0.74±0.02 


472 






Wavelength [ /um ] 



Fig. 4.— Spectrum of the star 2MASS J05351361-0518329 which is near dl32-1832. The 
model fit is shown as a gray line. No apparent water ice feature is detected in the spectrum. 
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Fig. 5. — Spatially resolved L-band spectra of d216-0939. Filled square shows K' photometric 
point, and gray lines are the best-fit model spectra to derive the water ice band for each 
positions. Right panels show the slit location and the area for extraction of the spectra, 
which are superposed on J, H, K' color composite image (upper) and L' image (lower) from 
the subtracted images in Figure O Gray box shows the slit location used on the first epoch 
of 2006-01-15-UT. 



Table 7. Spatially Resolved Water Ice Absorption in d2 16-0939 



Location 


Peak 






in d2 16-0939 


(/im) 




(cm^-*^) 


Northeast (NE) 


3.13 


0.45±0.02 


302 


Center 


3.10 


0.78±0.03 


479 


Southwest (SW) 


3.08 


0.57±0.04 


369 
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5. Discussion 

5.1. Water Ice in Individual Disks 

While d053-717 and dl41-1952 show no water ice feature at 3.1 /im band, moderate 
water ice absorption is detected for dl32-1832 and d216-0939. In the following, each 3/im 
spectrum is individually discussed. 



5.1.1. d053-717 



Although the morphology of th e silhouette disk f or this target obtained by HST indi- 
cates the nearly edge-on geomet ry (ISmith et al.l 120051) . the visual extinction has been esti- 
mated to be very low fAv=0.43: lHillenbrandl l 119971 ) from the optical spectroscopy. However, 
scattered light from central source may dominate the flux in the optical, which lead to the 
underestimate of the Ay. The object is located at moderate distance from O and B stars, 
and the UV environment is not too harsh to significantly affect the water ice distribution 
in the disk. Thus, the fact that no water ice feature at 3.1 /zm is detected toward d053-717 
is a direct evidence f or this object to have no significant amount of foreground extinction. 
Getman et al.l (120051 ) reported huge hydrogen column density {\og{NH)=22.7 cm~^) right 
at the position of d053-717 using the Chandra X-ray observatory and the X-ray source 
is likely embedded in the circumstellar disk. This column density is conve rted to the Ay 
of 31 .5 mag by using the relationship of Nh = (1.59 x 10^^) x Ay cm~^ (jimanishi et al. 



200ll ). which should produce significant optical depth of the water ice absorption at 3.1 /im. 
Therefore, it is concluded that the real source associated with the silhouette disk is located 
behind d053-717. In our AO images at K' and L', no other source than d053-717 is detected. 



5.1.2. dl 32- 1832 



The water ice absorption exists at 3.1 /im with an optical depth of 0.70 for this target. 
The equation between column density of water ice { NiH-jO)) and visual extinction (Ay) 
adopted in Paper 1, Ar(i720) = (1x10^^) x {Ay-Ay.o) Jwhittet et al.lll996l : lTeixeira fc Emerson 

19991 ) is used for estimate of Ay. First, the column density of the water ice toward dl32-1832 

is der ived to be 1.58 x 10^^ cm'^ using iV(if20) = ri^2o-Az//(2.0xl0-^^) fId'Hendecourt fc AUamandola 
19861 ). The threshold visual extinction of Ay.^Q ranges from 2 mag to 10 mag as an indicator of 
harshness of the UV field among various star forming clouds. As a result. Ay for dl32-1832 
is estimated to be 17.6-25.6 mag. This Ay is far greater than the average visual extinction 
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1-2 mag) of the foreground "neutral lid" in M42 (jO'Dell et al.lll992[ ) an d also significantly 
large r than the more local visual extinction (~10 mag) showed in Figure [T KjScandariato et al. 



201ll ). Therefore the detected water ice is concluded to be localized around the object dl32- 
1832. Since the nearby star 2MASS J05351361-0518329 has no water ice absorption, there 
is no common water ice material in front of those two stars. Assuming that dl32-1832 and 
2MASS J05351361-0518329 are located within the same region of M42, then the detected 
water ice originates from the region within a projected distance of 5'.'68 (2352 AU) around 
dl32-1832. 

Since the detailed analysis of the detected water ice profile for dl32-1832 is somewhat 
difficult with the S/N obtained, the overall proffie of the water ice feature extracted in 
Figure E] shows no broad absorption around 3.2 /xna and it is similar to the optical depth 
detected in the Taurus edge-on disks flTerada et al.l 120071 ). This suggests that the detected 
water ice grains originate from regions further out in the circumstellar disk rather than in the 
inner disk where large crystallized ice grains can be produced by grain growth and thermal 
processing. 



Interestingly, ISmith et al.l (120051 ) pointed out that dl32-1832 has been variable in the 
images taken with Ha filter of HST. The central star seems clearly visible right at the center 
of the silhouette disk on the first observing epoch of 1998-04-04-UT, and then it becomes 
invisible on the second epoch of 2004-01-18-UT. This variability suggests a change of the ex - 
tinction in the disk between two epochs or an Ha flare on the ffist epoch ( Smith et al. 12005 ). 



The i mage on another epoch (2005-04- 10-UT) provided as an online data by iRicci et al. 



(120081 ) still shows a little emission from the central star in the F658N and strong emission in 
the F775W and F850LP ffiters. The images taken with the F775W and F850LP ffiters show 
evidence of scattered light. It is unknown whether the central star was optically visible or 
invisible on the epoch of our observation (2011-10-23-UT). At the present time we cannot 
rule out variability due to an Ha flare and the amount of scattered light is not quantitatively 
understood. If the cause of the variability comes from a change in the extinction, then it 
would imply the water ice absorption would vary as well. 



5.1.3. dl41-1952 

While the silhouette disk around dl41-1952 is the smallest in size among our target 
samples, its disk mass derived from the submm measurement is the second most massive 
among the four targets (Table [1]). Assuming the abundance ratio between the water ice 
and the total dust mass to be similar among the protoplanetary disks, abundant water ice 
exists in the disk of dl41-1925. The silhouette disk shows the lowest inclination angle of 
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Fig. 6. — Normalized optical depth of the water ice feature for dl32-1832. In contrast to 
d216-0939, no 3.2 /xm enhancement of the optical depth is observed. 
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55-60°among the targets, and most likely our line of sight toward the central star does not 
pass through the circumstellar disk around the object. 



5.1.4. Spatially Resolved d2 16-0939 

As shown in Figure the d216-0939 spectra were spatially resolved with (T.'l resolution 
using the A0188 system with NGS mode. The bluish continua seen in the spectra of NE 
and SW positions are consistent with the bluer scattered light in those regions. The con- 
tamination from the central star mentioned in Section 4.2 affects NE and SW by the optical 
depth of 0.37 and 0.51, respectively. Therefore, the significant fraction of the optical depth 
measured at the locations of NE and SW originates from the water ice absorption at the 
center position, and the water ice absorption at NE and SW is very small (rjce=0.08 and 
0.06, respectively). These shallower optical depths at the position of NE and SW indicate 
that the water ice materials attributed to the absorption detected in NE and SW are located 
at greater height from the mid-plane. 

Focusing on the water ice feature detected in the center location of d216-0939, the 
profile of the water ice absorption is compared with the water ice profile reported in Paper 
I with a larger aperture (0'.'45 x O'.'S) in Figure [71 In the figure, the profiles of the water 
ice absorption are found to be almost identical with a possible decrease of the optical depth 
around 2.9 /im and 3.6 /im. The pure amorphous water ice has an absorption feature at 
2.9 /im , and therefore the decrease at 2.9 /im can be explained by the lower fraction of the 
amorphous water ice. Although the difference of the profiles at 2.9 yum is small, the water 
ice composition is quantitatively estimated in the following by assuming that the difference 
is real. We apply exactly the same method as used in Paper I, and the new result of 
the water ice optical depth in the region of 2.8-3.3 /xm is well matched to the composition 
of large-size crystallized water ice particles (77.2%±4.0%) and small-size amorphous water 
ice particles (22.8%±4.0%) with a minimum reduced and la uncertainty. Since the 
composition of the water ice detected on the first epoch is crystallized water ice particles 
of 71.4%±2.7% and amorphous water ice particles of 28.6%±2.7%, the results on those two 
epochs are overlapping within the la uncertainty. The optical depths of the crystallized 
water ice ijcTystai) are derived to be 0.55 and 0.69 on two observing epochs of 2006-01-15-UT 
and 2009-10-02-UT, respectively. Different apertures were used for spectroscopy on the two 
epochs, and therefore these derived optical depths cannot be directly compared. The large 
aperture of 0'.'45 x 0'.'8 used for extraction of the spectrum on the first epoch contains a large 
amount of scattered light from the central so urce. Since the large-size grains are thought to 



settle toward the mid-plane of the disk (e.g., IWilliams &: Ciezall201ll ). the spectrum of the 
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scattered light does not contain a significant amount of the large-size crystallized water ice. 
As a result, the scattered light component dilutes the optical depth of the crystallized water 
ice absorption. From the V image in Figure [H we estimate the fiuxes of the scattered light 
component (with an aperture of 0'.'45 x O'.'S) and the central region (with an aperture of 0'.'225 
X O'.'ll) of d216-0939 to be 2.29 x 10"^^ and 8.22 x lO'^^ W m'^ /im'^ respectively. If the 
scattered light is assumed to contain no large-size crystallized water ice feature, this dilution 
effect leads to the change of Tcrystai by 0.20. In reality, the scattered light includes at least 
partially the large-size crystallized water ice feature. Thus, the optical depth change due to 
the dilution effect is 0-0.20, and Tcrystai in the center region is estimated to be 0.55-0.75 on 
the first epoch. This optical depth range of the large-size crystallized water ice is consistent 
with Tcrystai =0.6 9 detected on the second epoch. 



5.2. No Variation of Large-size Crystallized Water Ice in d216-0939 Disk 



The variability of the water ice features has been detected toward sever al young stellar 
objects Jheiiiert et al.l l200ll: ISeck et al.lbood : ISecklboOTt iTerada et allboOTi ). In particular, 
Terada et al.l ( 120071 ) discovered the large variation (Axjce ~ 0.59) of the water ice in the edge- 
on disk around HV Tau C with a time interval of 2.32 years and derived the possible size of 
the water ice cloud to be 1.4 AU. No significant change of the large-size crystallized water 
ice feature in the time interval of 3.63 years indicates uniform distribution of the carrier of 
the water ice in the disk of d216-0939. Since the distribution of the crystallized water ice 
is not known well, we assume that it is located at the distance of 30 AU from the central 
star, which is suggested f or the first time as the lower limit of the crystallized water ice 
distribution in YLW 16A (ISchegerer fc WoliI 120101 ). Based on this assumption, the possible 
water ice cloud must be > 3.48 AU in size. 



5.3. Critical Parameters for Water Ice Detection 

Combined with the water ice features detected in Paper I, the relationship of the water 
ice optical depth with the inclination angle of the disks and the distance from the O and B 
stars is summarized in Figure [SI It is noted that the water ice for dl21-1925 is attributed to 
the foreground material. In addition, as pointed out in the above section, the lack of a water 
ice feature toward d053-717 is not well understood, but most likely d053-717 is not embedded 
in the silhouette disk. In this updated figure with the new samples, the rapid change of the 
optical depth can be found in the gray area for both the inclination angle of the disk and 
the distance from the O and B stars. Since UV radiation responsible for photodesorption 
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Fig. 7. — Comparison of the normalized optical depth of the water ice absorption between 
large aperture (0'.'45 x (/.'S) on 2006-01-15-UT and smaller region (0'.'225 x C'll) on 2009- 
10-02-UT around the photo center of d216-0939. 
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of the water ice is gradually changed as a function of the distance, it is difficult to explain 
such a sudden change only by the difference of the distance for dl41 -1925 (214^^18) and 



dl32 -1832 (294'.'19). On the other hand, as predicted in the model by (jPontoppidan et al. 



20051 ). the different inclination angle can produce this kind of rapid change of the water ice 
absorption. Therefore, we conclude the primary parameter for explaining the detected water 
ice absorption to be the inclination angle of the disk, and its critical angle is confirmed to 
be 65-75°. This suggests that the opening half-angles of the silhouette disks are around 
70°. This angle is a lmost the same angle derived for CRBR 2422-3423 in p Ophiuchus by 



Pontoppidan et al.l ( 120051 ). and the intense UV radiation does not significantly affect the 



critical inclination angle of the disks. 



6. Summary 

We present the AO images and spectra of the four silhouette disks in the outskirt region 
of M42 and M43, and the following results are reported. 

1. No water ice feature toward d053-717 and dl41-1953 was detected, and moderate 
absorption of the water ice was found in the disks of dl32-1832 and d216-0939. 

2. The water ice profile of dl32-1832 shows no signature of large-size crystallized water 
ice, and thus it may not have experienced grain growth in the disk. 

3. Water ice absorption in the d216-0939 disk was spatially resolved. The center region 
shows the deepest absorption, and the large-size crystallized water ice component is 
identical to that detected on 2006-01-15-UT. 

4. No time variation of the water ice in the d216-0939 disk was detected. Since the time 
interval is 3.63 years, this suggests a uniform distribution of the large-size crystallized 
water ice in the lateral direction of 3.48 AU. 

5. Summarizing the optical depths of the water ice absorption in the silhouette disks, no 
significant water ice absorption is detected for the disks with an inclination angle of < 
65°. Omitting the pecuhar sample of d053-717 (see Section 5.1.1), there seems to be a 
critical inclination angle of t he silhouette disks of 65-7 5° for water ice to be detected. 



which agrees to the model of iPontoppidan et al.l (120051 ). This implies that the opening 



half-angle of the silhouette disks is about 70°. 
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Fig. 8. — Relationship of the optical depth of the water ice absorption with the inclination 
angle of the disk (left panel) and the projected distance from O and B stars (right panel). 
Gray squares are from the result of Paper I, and black squares are reported here. The 
foreground water ice toward dl21-1925 is presented with open gray square. A sudden change 
of the optical depth is seen around the shaded area in both panels. 
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